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Abstract 

We extend our earlier macrostatistical treatment of hydro dynamical fluctuations about 
nonequilibrium steady states to viscous fluids. Since the scale dependence of the Navier- 
Stokes equations precludes the applicability of any infinite scale (hydrodynamical) limit, 
this has to based on the generic model of a large but finite system, rather than an infinite 
one. On this basis, together with the assumption of Onsager's regression hypothesis and 
conditions of local equilibrium and chaoticity, we show that the hydrodynamical fluctua- 
tions of a reservoir driven fluid about a nonequilibrium steady state execute a Gaussian 
Markov process that constitutes a mathematical structure for a generalised version of Lan- 
dau's fluctuating hydrodynamics and generically carries long range spatial correlations. 
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1. Introduction 



The theory of fluctuations is a key area of statistical physics, which is essential to both 
equilibrium [1] and nonequilibrium thermodynamics [2]. Further, Landau [3] introduced 
fluctuation theory into fluid dynamics by adding white noise terms to the Navier-Stokes 
equations and employing the fluctuation-dissipation theorem to relate their intensities to 
the viscosity and thermal conductivity coefficients. 

In relatively recent works [4, 5], we have presented a macrostatistical approach to 
the theory of hydrodynamical fluctuations of reservoir driven quantum systems about 
nonequilibrium steady states. This was designed to relate the stochastic properties of the 
fluctuations to the phenomenological continuum mechanics of these systems on the basis 
of general arguments centred on the hydrodynamical observables, subject to assumptions 
of macroscopic classicality, local equilibrium, chaoticity and a generalised version of On- 
sager's regression hypothesis* [2], all pertaining to states that may be far from global 
thermal equilibrium. On this basis, it was shown that, in a certain large scale limit, the 
fluctuations of the hydrodynamical variables executed a classical macroscopic stochastic 
process, whose parameters were expressed in terms of purely phenemenological quantities, 
namely thermodynamical variables and transport coefficients, the underlying quantum me- 
chanics of the system being buried in the forms of these variables. Among the results that 
ensued from a treatment of this stochastic process [4, 5] were a non-linear generalisation 
of Onsager's reciprocity relations and a proof that the spatial correlations of the hydro- 
dynamical variables are generically of long range. However, as it stands, the theory of 
those works is limited to systems whose phenomenological evolutions are scale invariant** 
and is therefore not applicable to viscous fluids, as described by the Navier-Stokes (NS) 
equations, since these are scale dependent. 

The object of the present article is to extend our macrostatistical treatment to viscous 
Navier-Stokes fluids and thereby to provide a mathematical structure for a generalised form 
of Landau's picture of hydrodynamical fluctuations and to establish that, generically, their 
spatial correlations in nonequilibrium steady states are of long range. However, the scale 
dependence of the NS equations precludes the applicability of a hydrodynamical limit***, 
which was basic to the methodology of the previous works. Consequently, the theory has 
to be based on the generic model of a large, but finite, system, rather than an infinite one. 

Our treatment here is based on a classical macrostatistical model, whose presumed re- 
lationship to an underlying quantum mechanics is only briefly indicated in the concluding 
Section. Thus, the model comprises a continuous distribution of matter that is confined to 

* In the present setting, this signifies essentially that small deviations of the hydrody- 
namical observables from their values in a nonequilibrium steady state evolve according to 
the same dynamical law whether this deviation arises from a spontaneous fluctuation or 
from a small external perturbation. 

** Specifically, those evolutions had to be invariant under space- time scale transforma- 
tions of the form x— >Xx, t— >X 2 t. 

*** This is a limit in which the space and time scales for the phenomenological description 
becomes infinite and the macroscopic equations of motion become exact. 



2 



a bounded spatial region Q and coupled to an array of reservoirs at its boundary. We as- 
sume that the scales of mass, distance and time whose units represent characteristic values 
of these variables are macroscopic: thus, for example, the constants % and ks of Planck and 
Boltzmann are extremely small on these scales*. We take the hydrodynamical observables 
to be the position and time dependent densities of energy, mass and momentum, and we 
assume that their phenomenological evolution is given by the NS equations. We then con- 
struct the model of the hydrodynamical fluctuations about a nonequilibrium steady state 
on the basis of assumptions of the Onsager regression hypothesis, local equilibrium and 
chaoticity. Although no hydrodynamical limit is available here, we exploit the fact that 
Boltzmann's constant is extremely small on the employed macroscopic scaling: specifically, 
we pass to a limit in which >0 in our formulation of local equilibrium conditions via 
Einstein's relation, P = const. exp(S , //cs) 5 between the equilibrium probability distribution 
of the macroscopic variables and the entropy function S. On this basis we obtain a gener- 
alisation of Landau's picture, wherein the hydrodynamical fluctuations execute a Gaussian 
Markov process whose parameters are completely determined by macroscopic variables in- 
volved in the phenomenological thermodynamic and hydrodynamic pictures of the system: 
the underlying quantum mechanics is assumed to be buried in the forms of these variables 
as functions of the control parameters. Futhermore, we show that the spatial correlations 
of the hydrodynamical variables are generically of long range. 

We present our treatment as follows. In Section 2 we formulate the thermodynamics 
and hydrodynamics of the model in purely phenomenological terms. In Section 3 we 
construct the stochastic process executed by the hydrodynamical fluctuations, subject to 
the assumptions specified in the previous paragraph. In Section 4 we prove that the spatial 
correlations of these fluctuations are generically of long range. We conclude in Section 5 
with some brief comments on the basis of the model and on its presumed relationship to 
its underlying quantum mechanics. There are two Appendices: the first is devoted to a 
calculation leading to a key formula, the second to a proof of a lemma. 

Note on distributions. As we shall represent the hydrodynamical fluctuations by dis- 
tributions, in the sense of L. Schwartz [6], we now specify our notations for these. We denote 
by V(Q), Vy{Q) and V T (Q) the Schwartz spaces of real valued, infinitely differentiable 
scalar, vector and second order tensor valued functions on the bounded open region O 
with support in that region. These spaces are reflexive and their duals are distributions, 
which we denote by V'(Q), V' V {Q) and V T (Q,), respectively. We define V{Q) and £>(fi) 
to be the Cartesian products V(Q) xV(Q) xV v (Q) and V V {Q) xV v (Q) xD T (0), respec- 
tively, and we denote their duals by V'{Vt) and V'(Q), repectively. Further, if a and ip 
are elements of a D-class space and its dual, respectively, then ip(a) will sometimes be 
denoted by (ip, a). We also employ angular brackets to denote the L 2 inner product (a, a') 
between pairs of elements a and ot of the same T>(Q) space. Evidently, these two uses of 
the angular brackets are mutually consistent since any V space is a subset of its dual. To 
avoid ambiguity in the definition of scalar products of tensor valued functions, we define 



* For example, if SI units and degrees Kelvin are appropriate, h and ks are of the order 
of 10~ 34 and 10~ 23 , respectively. 
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that of elements a (= {a.ij\i,j = (1, 2, . ., d)}) and a' (= {a^Ai^j = (1, 2, . ., d)}) of 2>r(0) 
to be 

(a, a') := / dxaij(x)a'^ i (x), (1.1) 

where, as elsewhere in this article, the repeated index summation convention has been 
employed. 

2. The Phenomenological Picture 

We assume that the fluid is a macroscopic system, E, that occupies a fixed bounded 
open connected region O of a <i-dimensional Euclidean space X, which constitutes the 
laboratory reference frame. We assume that E is in contact at its boundary, 90, with an 
array, 1Z, of reservoirs, and that these determine the boundary conditions for the flow 
of the system, in a way that we shall specify in Sec. 2.2. We further assume that its 
dynamics is Galilei covariant. We employ a continuum model for E, which we formulate 
on macroscopic scales wherein the magnitudes of its energy, mass and volume are of the 
order of unity: for simplicity, we take its volume to be unity. 

2.1. The Thermodynamic Potentials 

We assume that, at equilibrium, E is at rest, that its densities, e and p 0: of energy 
and mass are spatially uniform and that its entropy density is a function, sq, of these 
variables. Thus, as the volume of O is unity, eo, po and so( e o 7 po) are a ls° the total energy, 
mass and entropy, respectively, of E and satisfy the fundamental formula 

ds = f3(de - pdp ), (2.1) 

where is the temperature, in degrees Kelvin, say, not in units of ks, and \x is the 
chemical potential, as related to the mass*. Thus 

P=(P\ ( 2 - 2 ) 



and 



The pressure is then 



"=-^0'«» (2 - 3) 



p = /3 1 s (e ,p ) - e + /j,p . (2.4) 
and the heat function (enthalpy density) is 

e = e +p. (2.5) 



* We relate the chemical potential to mass, by Eq. (2.3), rather than particle number, 
since the continuum model does not involve any concept of the latter. 
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Turning now to the nonequilibrium situation, the densities of energy, mass and mo- 
mentum are generally non-uniform. We assume that they are locally conserved and we 
denote their local densities by e, p and j, respectively. These are functions of position 
x (GO) and time t (gR). The local drift velocity is defined to be 

u=j/p. (2.6) 

We denote the components of j and u, relative to some chosen coordinate system, by 
-Jd) and (ui,.., ltd), respectively. 

We define a local rest frame for the point x to be one that moves with velocity u{x) 
relative to the laboratory frame. It therefore follows from Galilei covariance that the local 
energy and mass densities relative to this frame are 

e = e - -pu 2 = e - — and p = p. (2.7) 
2 2p 

We assume that, even in a nonequilibrium state, the system is in local equilibrium in the 
phenomenological sense* that its local enthalpy density s is still given by Eq. (2.5) and 
its entropy density at a point x is the same function, so, of the energy and mass densities 
relative to a local rest frame as that governing the equilibrium entropy density relative to 
the laboratory frame. Thus, by the Galilei invariance of entropy [7 ], its density at the 
point x, relative to the laboratory frame, is the function s of e(x), p(x) and j(x) given by 
the formula 

s(e,p, j) = so(e ,Po), (2.8) 

We now compactify the notation by denoting the triple (e, p,j) by <fi. Thus, defining 
v := (d+2), <fi is the z/-component variable (0i, . ., 0„) := (e, p,ji, . ., jd) and s is a function 
of <f>. Its derivative s' (</>)(:= (ds/d<j>i, ..,ds/d^ v )) is then the conjugate, 6 (= (^i, .,#„)) 
of (p and, by Eqs. (2.2), (2.3), (2.7) and (2.8), it is given explicitly by the formula 

9 = = 0(1, -p + u 2 /2, - Ul , .., -u d ) (2.9) 

or, more compactly, 

6 = s'((j)) = (3(l,-p + u 2 /2,-u). (2.9a) 

The thermodynamic conjugate, tt, of s is the function of 9 defined by the formula 

n(9) = su^(s(<P)-9.<p). (2.10) 

We assume that all parts of this system are in a single thermodynamical phase at all times 
and consequently that the function s' is invertible and the supremum on the r.h.s. of Eq. 
(2.10) is attained when = [s']~ 1 (9). Hence 

7r(0) = s(0) - 94, with s'(0) = 9, (2.11) 



A further, macrostatistical kind of local equilibrium will be assumed in Section 3.3 
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which, together with Eqs. (2.2)-(2.4), (2.8) and (2.9a), implies that 



tt(0) = pp. (2.12) 

Furthermore, it follows from Eq. (2.11) that 

tt' (0) = -(j). (2.13) 

Consequently, the Hessians s"(<f>) (= (\d 2 s / d(f)jd(f)k\)) and it" (9) are related by the equa- 
tions 

n"{9)s"{4>) = s"(<P)n"(9) = -I, 

n"(9) = -s"(<P)- 1 . (2.14) 



2.2. The Navier-Stokes Equations. 

These hydrodynamical equations comprise local conservation laws for the energy, mass 
and momentum, together with constitutive equations for the associated fluxes. The local 
conservation laws are 

de dp dj 

— + V.q = 0; -£ + V.j = 0; and -± + V.r = 0, (2.15) 

where q is the energy current, r is the stress tensor [tm] and (V.r)/- := drki/dxi. The 
constitutive equations are then [3] 

q= {e + ]-pu 2 )u- o.u + kV(3, (2.16) 

t = pi + puu — a (2-17) 

and 

a = 7i (Du - 2d-\V.u)l) + 7 2 (V.u)J, (2.18) 

where 

fn \ du k . dui 

{Du)kl = ^ + toT k > (2 ' 19) 

uu is the dyadic whose kl component is UkUi and k, 71 and 72 are positive, scalar valued 
functions of (3 and p that represent the thermal conductivity, the bulk viscosity and the 
shear viscosity, respectively. 

Thus, in view of the one-to-one correspondence between the variables (e,p,j) and 
(f3,p,,u), the combination of the conservation laws (2.15) and the constitutive equations 
(2.16)-(2.18) describe an autonomous evolution of the hydrodynamical variables. Further, 
as remarked at the beginning of Section 2, the boundary conditions are determined by the 
reservoirs with which the system is in contact. Specifically, we assume that the contact 
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is restricted to the surface, 90, of O and the boundary conditions are taken to be the 
following ones. 

(a) u = on 90. 

(b) At any point of 90 that is in contact with a reservoir, the values of the temperature 
and chemical potential of E are just those of the reservoir. 

(c) Any part of 90 not in contact with a reservoir is insulated and the normal part of V6> 
is zero there. 

We now compactify the hydro dynamical equations by expressing them in terms of 
the triples <fi '■= ( e ?P?j) an d X '■= ( c Li3i T )i the latter of which we shall term the flux. 
Inserting subscripts t to indicate time dependence of the hydro dynamical variables, the 
local conservation laws (2.15) may be expressed as the single equation 

^ + V. X t = 0, (2.20) 

and Eqs. (2.16)-(2.18) comprise a constitutive equation for xt of the form 

Xt (x) = G{(j>t : x). (2.21) 
By Eqs. (2.20) and (2.21), (j) t evolves according to the autonomous equation of motion 

d<p t (x) 



dt 

where 



= ^(0 t ;x), (2.22) 



T{(j>ux) :=-V.g{(j>ux). (2.23) 

We assume that the equation of motion (2.22) has a unique stationary solution, subject to 
the prevailing boundary conditions, described above. We denote this solution by <j)(x) = 
(e(x) : p(x) : j(x) and denote the corresponding steady state value of 9(x) by 9{x). 

The Perturbed flow. We assume that the stationary flow is stable under 'small' pertur- 
bations 5(pt(x) = (5et(x), 5pt(x), 5jt(x)) of </>(rr), that preserve the boundary conditions. 
The resultant linearised equation of motion is then 



-Q t Ht{x) = (C5cp t )(x) := — F@ + \6cf> t ;x) lx=0 . (2.24) 



Thus, by Eqs. (2.23) and (2.24), 

C5<p t = -VXdfa, (2.25) 

where 

7 



= ^rG(<P + X5<p t ;x), x=0 . (2.26) 



Further, by Eqs. (2.21) and (2.26), the increment 5\t in Xt is given by the formula 



5 X t(x)=>C5fa(x). (2-27) 

In anticipation of the demands of the macrostatistical picture of the system, formu- 
lated in the next Section, we assume that 5 fa is a distribution, in the sense of L. Schwartz 
[6]. Thus, since 8 fa = (Se t ,8p t ,Sj t ), where the first two components are scalar fields and 
the third is a vector field, we assume that 8 fa is an element of the space X>'(0), defined at 
the end of Section 1. Thus, employing the notations introduced there, it follows from the 
definition of 8 fa that, for F = (f,g, h), with f,g&V(Q) and fteXV(f2), 

Sfa(F) = Set(f) + Sp t (g) + Sj t (h), (2.28) 

where the three terms on the r.h.s. are the integrals of Set, 8pt and Sjt against /, g and h, 
respectively. 

We now assume that the linear operator £, defined in Eq. (2.24), is the generator of 
a one parameter semigroup T := {T t \teH + } of linear transformations of V'(Q). Hence, by 
Eq. (2.24), 

Sfa = T t 5(po V tGR + , (2.29) 

from which it follows that 

8fa(F) = S(j)(T*F) V FeV(Q). teK+, (2.30) 

where T* := {T/|tGR+} is the dual of the semigroup T: its generator, £*, is the dual 
of C We assume the dissipativity condition that, for all perturbations 5(f), T t S(j) tends to 
zero, in the V'{Vt) topology, as t— >oo. As the space V{Vt) is reflexive, this is equivalent to 
the assumption that 

V(fl) : lim^ooT/F = V FeP(fi). (2.31) 

Re-expression of NS in terms of 9. Since the constitutive equations (2.16)-(2.19) 
relate the heat and mass currents, as well as the stress tensor, directly to 9, it is sometimes 
convenient to express the r.h.s. 's of Eqs. (2.22)-(2.24) in terms of this variable rather than 
<j>. Thus we rewrite Eq. (2.22) as 

^ = -V.*(0 t ;s), (2.32) 
where ^{9 t ;x) = T{fa;x). Eq. (2.24) for the perturbed flow then becomes 

-Sfa=AS9 t := -—v.*(6 + \60f,x)\x=o. (2.33) 

Further, by Eq. (2.13) and the equivalence of Eqs. (2.24) and (2.33), 

A59 = C8(f> = -Cn"(9)59 
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for all perturbations 66 of 6 and hence 

A = -Cn"(6). (2.34) 

It follows immediately from this equation and the symmetry of rc(0)" that the dual, A*, 
of A is related to that, £*, of £ by the formula 

A* = -7r"(0)£*. (2.35) 

In particular, as we shall show in Appendix A, it follows from our definitions that the 
equilibrium form, A* q , of A* is given by the equation 

K q (f,9,h) = 

-/r^v/ - /rVv<7 + /S-SiV.cd/i - 2d-\v.h)i) + /rSvp?./*)) 

V(f,g,h)eV)Q). (2.36) 

3. The Stochastic Fluctuation Process 

According to elementary statistical mechanics, the thermodynamic and hydrodynamic 
variables undergo fluctuations, which are not taken into account in the phenomenological 
picture of the previous Section. We now seek to provide a general description of the 
hydro dynamical fluctuations about nonequilibrium steady states by treating the fields 
4> = (e,p,j) and x = {o.iji r ) as expectation values of random fields <fi = (e,p,j) and 
X = (QiJi t), respectively, where the tensor f , like the phenomenological r, is symmetric*. 
We assume that these fields satisfy the local conservation law given by the canonical 
analogue of Eq. (2.20), namely 

^ + V.xt = 0. (3.1) 

The differences between the random fields <f>t an d Xt and their classical expectation values 
then represent the hydro dynamical fluctuations. In a standard way, we normalise them by 
a factor A/" 1 / 2 , where M is chosen to be the ratio of characteristic values of corresponding 
macroscopic and microscopic quantities. In the present situation, where we are formulat- 
ing the model on a macroscopic scale, a natural choice for J\f is the reciprocal, fc^ 1 , of 
Boltzmann's constant, which arises in Einstein's formula, P = const. exp(5 , //cs), relating 
the equilibrium probability distribution of the macroscopic variables to the entropy S. 
Thus, we define the fields representing the fluctuations of <f> t and xt about their steady 
state values to be 

£ t (x) = k- 1/2 (Mx)-0) (3.2) 



* This symmetry property may be regarded vas basic, as it prevails in standard micr- 
roscopic pictures of stress tensors. 
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and 

rit(x) = k^ 2 (xt(x)-x), (3-3) 

respectively. Hence, denning 

Ct, s ■= / dur] u , (3.4) 

J s 

it follows from Eqs. (3.1)-(3.4) that 

6 - is = -V.Ct )fl V t, seK. (3.5) 

In accordance with the general requirements of field theories [8] , we assume that these 
fields are distributions, in the sense of L. Schwartz [6]. Specifically, in the notation of the 
last Note of Section 1, we assume that £ t e£>'(0) and r} t &V'(Q). We denote the smeared 
fields obtained by integrating £ t and £t,s against test functions F (G£>(fi)) and G (eX>(0)) 
by £t(F) and rjt(G), respectively; 

We aim to derive the stochastic process executed by the fluctuation fields £ and £ 
from the assumptions of 

(a) a generalised version of Onsager's regression hypothesis, 

(b) a macrostatistical local equilibrium hypothesis and 

(c) a chaoticty hypothesis for the random currents and stresses. 

Further, as the continuum model harbours only macroscopic variables, we take any coun- 
terparts to microscopic correlation lengths and memory times of this model to be zero. 

3.1. Regression Hypothesis. 

The regression hypothesis signifies that the fluctuations of the hydrodynamical ob- 
servables evolve, from a given starting point, according to the same dynamical law that 
governs the evolution of small perturbations of those variables from their steady state 
values. Thus, noting the formula (2.30) for the latter evolution and denoting conditional 
expectations of the stochastic variables with respect to the value of the field £ at time t 
by E(.\£ to )i we assume that 

E(Z t (F)\Z t0 ) = Z t0 ((T*(t-t )F)VFeT>(n), t ,t(>t )eR. (3.6) 

Since the process £t is stationary, an immediate consequence of this formula is that 

E(Zt(F)b(F')) =E(t[T*(t-t')F}t(F')) VF,F'eP(ft), t,t'(<t)eR. (3.7) 

We now note that, by Eq. (2.27), the increment in the integrated phenomenologi- 
cal flux, duxu, due to a perturbation 5(f) of <fi is J* dulC5<j) u . Correspondingly, as the 
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phenomeno logical dynamics of the model is secular, we designate the secular part of the 
integrated fluctuation flux, £ tjS , to be 

fduKi u . (3.8) 

J s 

We then define the remaining part, £ tjS , of Ct, s to be the stochastic part of the integrated 
flux. Thus 

Ct,a := Cm - / du)C£ u . (3.9) 



We shall presently show, in Prop. 3.1, that this field does indeed enjoy strong, Wiener-like, 
stochastic properties. 

By Eqs. (2.25), (3.5) and (3.9), 

£ t (F)-£ s (F)= I duU£*F) + w t , s (F) V s,teR, FeV(Q), (3.10) 

J s 



l.C. 



£t - 6 = / duCi u + w t , a , (3.10a) 

J s 

where 

w M (F) = C M (VF). (3.11) 

The following Proposition, which was proved in [5, Prop. 5.1], establishes that w simulates 
a Wiener process, at least as far as its two-point function is concerned; and hence that Eq. 
(3.10a) is an integrated Langevin equation. 

Proposition 3.1. Under the assumption of the regression hypothesis, w has the 
following properties. 

(i) 

E(w tiS (F)£ u (F')) = V t>s>u, F, F'eV(Q). (3.12) 

(ii) 

V s,t(>s), s',t'(>s) GR, F,F'eV(tt), (3.13) 
where \I\ denotes the length of an interval I in R. 



3.2. The Chaoticity Hypothesis. 

We now strengthen Prop. 3.1 by the assumption that the stochastic part, C of the 
integrated fluctuation flux is chaotic in the sense that its space-time correlations are of 
microscopic range, idealised here as zero range on our macroscopic scale. This assumption 
is designed to represent Boltzmann's molecular chaos hypothesis, as transferred to the 
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stochastic flux. Here we take it to signify that the field ( is Gaussian, due to statistical 
independence of its values in disjoint space-time regions. Thus, our chaoticity hypothesis 
comprises the following conditions. 

(C.l) The process ( is Gaussian, and 

(C.2) E (( Sjt (G)( s / jt > (G')) vanishes if the intersection of either the intervals [s, t] and [s', t'] 
or of the supports of G and G' is empty. 

Further, we supplement these conditions by the following continuity assumption. 

C. The correlation function of (C.2) is continuous in its time variables. 

The following Proposition was proved in [5, Prop. 5.2]. 

Proposition 3.2. Under the assumption of (G2) andC, the two=point function of 
( takes the following form. 

E(Ct, s (G)Cs',t'(G')) =T(G,G')\[s,t)n[s',t'}\ VG,G'eP'(fi), s^s'JelL, (3.14) 

where V is a continuous bilinear form on Z>(0)(g)Z>(0), whose support lies in the region 
{(x, x')eVt 2 \x = x'}. 

We shall derive the explicit form of V from the local equilibrium hypothesis in Section 
3.3. The following Corollary to Prop. 3.2 is a simple consequence of Eq. (3.13), assumption 
(C.l) and Prop. (3.2). 

Corollary 3.3. Under the above assumptions, the process w is Gaussian with zero 
mean and its two-point function takes the form 

E(w tyS (F)w t ,, s ,(F')) =T(VF,VF')\[s,t}n[s , ,t'}\ VF,F'eP(Q), t,s,t',s'eK. (3.15) 



The following Proposition was shown in [5, Prop. 5.5] to ensue from the Langevin 
equation (3.10a), Prop. 3.1 (i) and Cor. 3.3. 

Proposition 3.4. Under the above assumptions, £ is a Gaussian Markov process, 
and the fields w tjS and £, u are statistically independent of one another if s and t are greater 
than or equal to u. 

3.3 Equilibrium and Local Equilibrium Conditions. 

Equilibrium Statistics o/£. We base our formulation of these statistics on the canonical 
version of Einstein's formula for the probability distribution, P of the hydro dynamical 
variables. This is given formally by the equation 

P e ,=co nS t.expfe> / -**(*)]), 

JO. 
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i.e., by Eq. (3,2), 



p 

1 eq 



= const.exp^ 1 I dx[s(</> + k^ 2 ^{x)) - 9.(</> + k^ 2 ^{x)]). (3.16) 



Since the integrand of this formula is maximised at £ = 0, its ratio to ks reduces to 
(£(x).s"(@)£(x))/2 in the limit k B ^0, which we term the Botzmann limit: here the dot 
denotes the scalar product. Thus, operating henceforth in this limit, the equilibrium 
characteristic function for the fluctuation field £ is given formally by 

F eq (exp [i£ (F )}) = const. J ^ [ i Uil ; , V ; V _ , ' V IJi V FeV(Q), 

(3-17) 

where X>£ denotes functional integration w.r.t. the field £. This formula may properly be 
defined by resolving O into a set of cells, A J, denoting the values of £ and F at the centre 
of AJ by £j and Fj, respectively, and then expressing Eq. (3.17) as 

F eq exp[^(F)] = hm A nj - v '- , (3.17a) 

where limA is the limit in which the cells shrink to points. It now follows easily from 
Eqs.(2.14) and (3.17a) that, in the notation specified at the end of Section 1, 

£eq(exp[i£(F)]) = exp(-^(F,7r"(#)F» V FeV(Q). (3.18) 

Thus, at equilibrium, £ is a Gaussian random field, with zero mean and two- point function 
given by the formula 

E eq {£(F)) 2 ) = 1 -(F,n"(e)F)V FeV(n), 

or equivalently, by polarisation, 

F eq (£(F)£(F')) = (F, 7r"(0)F')> V F, F'eP(Q). (3.19) 



Equilibrium two-point function for ( . By Eqs. (2.34), (3.11), (3.13) and (3.19), 
F eq (C M (VF)C^(F')) = [(A: q F,F') + (F,A: q F')]|[ S ,t]n[ S ',t']| 

V F,F'eV(Q), s,t,s',t'eK. (3.20) 
Hence, for F = (/, o, h) and F' = (f',g',ti), it follows from Eqs. (2.36) and (3.20) that 

F eq (C M (VF)C f ', s '(VF0) = 
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[2k(V/, V/') + /3-Si " 2eT 1 (V./i, V./i')) + 2/3- 1 72 (V./i, V./i')] x 

|[5,*]D[s',*']|. (3.21) 

We now express £ t>a in the form 

Ct,s = (CtJ i CtJ i CtJ)i (3.22) 

where the components are the time integrals of the stochastic fluctuations of q,j, r and 
therefore lie in T>' V (Q), V' V (Q) and T>' T (Q), respectively. Thus, 

Cm(vf) = cg(v/) + cS(v^) + c t (3 s } (v/i) 

and similarly 

C m (vf') = cg/W) + cff (W) + cffCW). 

The substitution of these last two equations in Eq. (3.21) yields the formulae 

MS^/K'i'W')) = 2«(v/, v/')|[M]n[ a ',f]|, (3.23) 

and 

E eq {C^(Vh)C^(Vh')) = [2r 1 li{(Dh,Dti)-2d- 1 (V.KV.ti))+r 1 l2{V.KV.ti))] 

\[s,t]n[s',t'}\. (3.24) 

and implies that all other two-point functions for the components of (^ jS (V-F) are zero 

The following lemma for the generalised functions and £^ on O corresponding 
to the distributions denoted by the same symbols will proved in Appendix B. 

Lemma 3.5. Under the above assumptions, supplemented by the condition that the 
equilibrium two-point function of ( is locally translationally and rotationally invariant, and 
indicating the components of X -vectors by subscripts i,j,l,m 

M3? ; «(*)$lW*0) = 

2/3 -1 (7i(<%^m + S im S jm - 2d~ 1 5 i i5 jm 5 j i) + i28u8 jm )8{x - x')|[s, t]n[s' \t']\ (3.26) 

and all other two-point functions of the components of ( t ,s o,re zero. This result concurs 
with that of Landau [3, Eqs. (132.11-13)]. 

In order to re-express the two-point functions of ( in terms of smeared fields, we denote 
elements G and G' of X>(fi) by triples (a, 6, c) and (a',b',c'), respectively, where the first 
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two components of both G and G' lie in XV (O) and the third lie in XV (fi). We then define 
the elements c /(1) of X> T (0) and c (2) , c'^ of X>(fi) by the formulae 

:= Cij + Cji] c£j ■* := + (3.27) 

and 

„(2) ._„... „/(2) _ , 

"J J' 

and infer from Eqs. (3.26)-(3.28) that 



c^:= CjJ ; c'^=c' (3.28) 



E eq (Ct,s(G)Ct',s'(G , )= [2 K (a,a')+ 

/3- 1 7i(c (1) ,c'W) + 2/3- 1 (72-2rf- 1 7l )(c (2) ,c'( 2 ))]|[ S ,t]n[ S / ,t / ]|. (3.29) 

This formula and Eq. (3.19) constitute our equilibrium conditions. In order to obtain 
their local properties, we consider their forms when their test functions are concentrated 
around points of O. Thus, for xo<EQ and e€R+, we define the transformations F—^F XOiC 
and G— >G XOi£ of X>(fi) and X>(fi), respectively, by the equations 



x ,e 



e - d / 2 F(e- 1 (x-x )) (3.30) 
and 

G x ^ = e- d ' 2 G{e-\x-x )) (3.31). 

We then remark that Eqs. (3.19) and (3.29) are invariant under these transformations and 
therefore that they enjoy the local (punctual!) property that 

lim^eq^^xo,,)^^,,)) = (F , n"(6)F'} V xoGO, F , F'eV(Q). (3.32) 

and 

lim^o^eq (Ct,s(Gx ,e)Ct',s' (G£ , e )) = 2 K {a,a') + 

/3- 1 7i(c (1) ,c / «) + 2/3- 1 (72-2rf- 1 7l )(c( 2 ),c^)]|[ S ,t]n[ S ',t']|. (3.33) 

Local Equilibrium Conditions. We now assume that, even in a nonequilibrium steady 
state, the two point functions of £ and ( enjoy the same local properties as at equilibrium. 
Thus, bearing in mind that (3, 9, ft, 71 and 72, are generally position dependent in the 
nonequilibrium situation, we take the local equilibrium conditions to be the following ones. 

\im ei0 E(^F Xo ^(F XOj€ )) = (F , n"(6(x ))F') V x GO, F , F'eV(Q). (3.34) 

and 

\im eio E(C t , a {G XOie )Ct',s'{G' XOie )) = [2 K (x )(a,a'} + 

/3(*orS(*o)(c ( V (1) >+2/3(£o)-^^^^ (3.35) 
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The following Proposition, whose proof ensues from a trivial modification of that of 
[5, Prop. 5.3], provides an explicit formula for the bilinear form T, which governs the form 
of the two-point function for ( according to Eq. (3.14). 

Proposition 3.6. Under the assumptions of Prop. 3.2, together with the local equi- 
librium condition (3.35), Y is given by the following formula. 

T(G, G') = 2{a,Ka') + (c«, /T^ic'^) + 2(c^ 2 \ /?" 1 ( 72 - 2cT 1 7l )c /(2) >, (3.36) 

where a, a', c^, c'^ c*- 2 -*, c'^ are related to G and G' according to the above specifications 
and now (3 1 , k, 71 are functions of position, through their dependence on f3 and fx, that 
act multiplicatively on V(Q). 

3.4.The Macrostatistical Model. 

The stationary processes £ and w, which are connected by the integrated Langevin 
equation (3.10a), comprise our macrostatistical model. In view of Cor. (3.3) and Prop. 
(3.4), these processes are both Gaussian, and the two-point function of w is given by Eqs. 
(3.15) and (3.36). To complete the formulation of the model, it remains for us to obtain 
the two-point function of the process £. 

To this end, we infer from Eq. (3.10a) that, since C is the generator of the semigroup 

T, 

& = T t _ t( Xt + [ T t _ s dw sM V t>t 

J to 

and hence that 

6(F) = 6 (T t *_ to F) + f dw s , t0 (T t *_ s F) V FeV(n), t,t (<t)eR. (3.37) 

In view of the stationarity of the ^-process, it follows from Eq. (3.37) and Prop. 3.4 that 
the static two-point function for the field £ is 

W(F,F') := E(t(F)Z(F') = E^(T t % o F)^T t % o F') + 
[ dsT(VT*F,VT*F')V F,F'eV(Q), t,t {<t)eR (3.38) 

J t 

On invoking the dissipative condition (2.31) and passing to the limiting form of Eq. (3.38) 
as to - ^ — 00, we obtain the formula 

POO 

W(F,F')= dsT(VT*F, VT*F') V F, F'ef>(0,). (3.39) 
Jo 

We note that, as C* is the generator of the semigroup T*, it follows from this formula and 
Eq. (3.38) that 

W(£*F, F') + W(F, C*F') + r(VF, VF') = 0. (3.40) 
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Moreover, in view of the stationarity of the process £, it follows from Eqs. (3.38) and 
(3.39) that the two-point space time correlation function is given by the formula 

/•OO 

E(Zt(F)&(F') = / dsr(VT*_ tl+a F, VTZF') V F,F'ef>(Q), t,t'(<t)eR. (3.41) 
Jo 

As the process £ is Gaussian, this completes our formulation of the model. 

4. Long Range Correlations. 

We term the space correlations of £ to be of short range, which we idealise as zero 
range in our macroscopic scaling, if the support of its static two-point function W lies in 
the region {(x, x')<EQ 2 \x = x'}. Then, by direct analogy with the proof of [5, Cor. 5.4], it 
follows from Schwartz's point and compact support theorems [6, Ths. 35 and 26] that if 
the correlations are of short range then W takes the form 

W(F, F') = (F, n"(9)F') V F, F'eV(Q). (4.1) 

On the other hand, we term the static space correlations of £ to be of long range if this 
condition is violated. Thus, in the present context, 'long' is taken to mean non-zero on 
the employed macroscopic scale. 

Our aim now is to establish that these correlations are generically of long range, i.e. 
that it is only in exceptional circumstances that the condition (4.1) is valid. To this end, 
we note that, in view of Eqs. (2.35) and (3.40), this condition may be expresses in the 
form 

(F, A*F') + (F\ A*F) = r(VF, VF') (4.2) 

. Thus, the condition for long range correlations is that of the violation of Eq. (4.2) for 
some F and F' in X>(fi). 

Proposition 4.1. A sufficient condition for the process £ to have long range space 
correlations is that one of the following ones is violated. 

u = (4.3) 

and 

v.((/rV«A.-«0W) = o, (4.4) 

where Kp and are the derivatives of k w.r.t. (3 and fx, respectively. Since these conditions 
can be satisfied only by certain particular forms of the space- dependent variable (3, fi and 
u, this signifies that the space correlations of £ are generically of long range. 

Remark. It will be seen that the proof of this Proposition is based on the choice 
F = F' = (/, 0, 0) for the forms of the test functions, with the result that long range 
correlations prevail if either the condition (4.3) for the drift velocity or (4.4) for the thermal 
conductivity is violated. We remark here that further sufficient conditions for long range 
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correlations, expressed in terms of the bulk and shear viscosities, may similarly be derived 
from other choices of the forms of F and F' . 

Proof of Prop. 4.1. It suffices to show that the condition (4.2) for zero range spatial 
correlations implies Eqs. (4.3) and (4.4). To this end, we choose both F and F' to be 
(/, 0, 0) and infer from Eq. (3.36) that 

T(VF,VF) = 2 I dxK(Vf) 2 . (4.5) 
Jo. 

On the other hand, by Eqs. (2.32) and (2.33), 

(F, A*F') = (AF, F') = ^(e + XF),VF') lx=0 V F,F'<=V(Sl). (4.6) 

In order to treat this formula for the case where F = F' = (/, 0, 0), we define 

Xf (x) :=e + X(f(x),0,0) (4.7) 

On asserting the ^-dependence of the energy current q by referring to q(x) as q(9;x) — , 
defining 

q Xf (x) :=q{0 X f,x) (4.8) 

and noting that, by identification of Eq. (2.15) with Eq. (2.32), \I/ := (q,j,r), we infer 
from Eq. (4.6) that 

(F,A*F) = [-^^dxVf(x).q Xf (x)} ]x=0 . (4.9) 

In order to determine the explicit form of q\f in terms of the position dependent variables 
(3, \i and u, we define their canonical counterparts (3\f, nxf and u\f, respectively, that 
correspond to 9\f by the version of Eq. (2.9) obtained by imposing the subscript Xf to 
each of its terms. Thus 

Oxf = /5a/ (1, -fixf + \u\ p -u Xf ). (4.10) 
It then follows from Eqs. (2.9) and (4.10) that 

Pxf = P + A/, iixf = (1 + A/?" 1 /)" V - A/3/(l + A/3- 1 /)- 2 ^ 2 and 

ux f = (l + \r 1 f)~ 1 u. (4.11) 

Further, since the energy current, q, defined by Eq. (2.16), is a functional of 9, we express 
qxf as the corresponding functional of 6\f, i.e. of /3a/, A*a/ an d u\f. Thus, 

Q\f = e(P\ f , n\f)u Xf + -p(Pxf, LixfWx f ux- 

7i (/^A/, n\f) {Du X f - 2d~ 1 V.uxfl).u X f - 72 (/3a/, //a/)(V.«a/)«a/ + ^(P\f, ^a/)V/3 A /. 

(4.12) 
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Since, by Eq. (4.11), 

(il%=" = /; (^f),.=o = -/'-'(" +« 2 )/; = (4.13) 

it follows from Eqs. (4.11)-(4.13) that 

(^f ) |a=o = + 0~ V ( uV f + - Zd-^.C/Ti)/)^ + ( 72 («.V/)«+ 

( Kf) - (3~ V«„)(V/?)/ + ^; m)/, (4.14) 

where and are the derivatives of k w.r.t. (3 and /z, respectively, and C is a vector 
valued functional of (3, \i and u that vanishes when u = 0. 

It follows now from Eqs. (4.5) and (4.14) that the condition (4.2) for short range 
correlations reduces to the formula 



/ 



ete/T 1 [7i ((W) V + (1 - 2d- 1 )(u. V/) 2 ) + 72 (u. V/) 2 )] + 

/ dx^/ 2 V.((/3-V«/i " «/?)V/3 - C) = V /eP(fi). (4.15) 
In particular, if we replace / here by f Xo ,e where xo&Q, eGR+ and 

fx ,e = €^-^ 2 f(e- 1 (x-x )), (4.16) 
then the passage to the limit e— >0 annihilates the second integral and yields the equation 

dxPixo)- 1 [ 7 i(xo)((V/(x)) 2 W (x ) 2 + (1 - 2d- 1 )(u(x ).Vf(x)) 2 ) + 



/ 

Jo. 



72(^0) (u(x ).Vf(x)) 2 } = V x GO, feV(fl). (4.17) 

For d>2, this implies that the velocity field u vanishes and consequently that the condition 
(4.15) reduces to that of the vanishing of the second integral of that equation, with C = 0; 
and for d = 1, 71 may be equated to zero in the NS equation (2.18) and therefore the 
same conclusion is valid. Thus, the short range correlation condition implies Eqs. (4.3) 
and (4.4), and therefore the violation of either of those formulae implies that the £ process 
has long range spatial correlations 



5. Concluding Remarks 

We have shown that the stochastic process executed by the hydrodynamical fluctua- 
tions of the continuum model about a nonequilibrium steady state is completely determined 
by the conditions of Onsagerian regression, local stability and chaoticity. This process, fully 
specified in Section 3.4, constitutes a mathematical generalisation of Landau's picture of 
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hydrodynamical fluctuations. In particular, by Prop. 4.1, the process generically carries 
long range spatial correlations. This result appears to be new within the framework of the 
Navier-Stokes equations, though it has been previously suggested on rather general heuris- 
tic grounds [9, 10] and proved for a certain classical stochastic (non-Hamiltonian) model 
[11]-[13]. Most notably, it marks an important difference between equilibrium and nonequi- 
librium properties of hydrodynamical fluctuations, as the spatial correlations of the former 
are of short (microscopic) range, except at critical points. We remark here that there is no 
corresponding qualitative difference between the time correlations of the hydrodynamical 
fluctuations about equilibrium and nonequilibrium, since the regression hypothesis implies 
that the time scales of both are the macroscopic ones of the Navier-Stokes flow. 

The treatment of hydrodynamical fluctuations in this article has been based on a 
classical macroscopic continuum model, E, of a fluid. Presumably this should arise from 
an underlying quantum mechanics, at the microscopic level, in the following way. One 
assumes that the quantum system, E qu , consists of N particles of one species that is 
confined to the region O and coupled at the boundary <9fi to an array, 1Z, of reservoirs 
whose temperatures and chemical potentials are just those of the macrostatistical model. 
Assuming that these control variables are not the same for all the reservoirs, the system 
(S qu + 71) will evolve, under rather general conditions, to a nonequilibrium steady state u 
[14]. The hydrodynamical observables of E qu may then be formulated along the lines of [4, 
5] in terms of the natural counterparts <^ qu and x qu , of the random classical fields (p and 
X of the model E, though those quantum fields are now operator valued functions of the 
positions and momenta of the particles of the system. We denote their evolutes at time £, 
as governed by the dynamics of the composite (E qu + 7?.) by </> qu; t and x qu; t, and we assume 
that these are distributions of class V(Q) and T>'(Q) and that their expectation values for 
the state ui are the classical steady state fields <p and X, respectively, of Section 2. We 
then define the quantum fluctuation fields ^ qu -t and Cqu;t,s by the canonical analogues of 
Eqs. (3.2) and (3.4) and denote the smeared fields obtained by integrating them against 
test functions F (eP(fi)) and G (eP(fi)) by £ qu; t(-F) and Cqu;t, a (<3), respectively. The 
correlation functions given by the expectation values, for the state u of the monomials in 
the £ qu; t(-F)'s and Cqu;t,s(G r )'s then represent the quantum stochastic process [15] executed 
by the hydrodynamical fluctuations. Further, under a condition of macroscopic classicality, 
whereby the correlation functions are invariant under reordering of the constituent smeared 
fields, this process is classical. Thus, under the assumption that this condition is fulfilled, 
possibly up to corrections that are o(l) w.r.t. h, ks and microscopic relaxation times 
and attenuation lengths, the fluctuation process simulates a classical one The further 
assumptions of Onsagerian regression, local stability and chaoticity then lead precisely 
to the classical macrostatistical one presented here. The ultimate test of the physical 
validity of that model is that its correlation functions are those of the hydrodynamical 
observables of E qu , up to the above microscopic corrections 

Appendix A: Derivation of the Formula (2.36) 

We assume that, at equilibrium, u = and (3 and \i are spatially uniform. The same 
is therefore true of e, p, k, 71 and 72, since these are functions of the latter two variables. 
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Thus, by Eq. (2.9a), 

0=(/?,-AM) (Al) 

and, as (p = (e, p, pu), it follows from Eq. (2.33), together with the identification of the 
NS equations (2.15)-(2.19) with Eq. (2.32), that the equilibrium form, A eq , of A is given 
by the formula 

A eq 5(/3, -/3/J,, -f3u) = 

(sV.Su - nASfi, -pV.Su, -ppV5(3 - p^Vdp + 71 V. (Ddu - 2d~ 1 (V.du)I) + 7 2 V(V.fo/)) , 

(A2) 

where pp and p M are the derivatives of p w.r.t. j3 and p, respectively. Since 4> = (e, p, pu), 
it follows from Eqs. (2.5) and (2.10)-(2.12) that 

Pp = (e - pp) and p^ = p. (A3) 

Hence, by Eqs. (A. 2) and (A. 3), 

K eq 50=K eq (5l3, -/35p - p,8/3, -/35u) = 

(-eV.Su - kA<5/3, -pW.Su, /0 -1 (e - pp)V5(3 - pVSp+ 

7iV.(D5w - 2d- 1 (V.Su)I) + 7 2 V(V.<5u)). (A4) 
Equivalents, defining (t^ 1 ), 8d^ 2 \ 56^) := 56 = (5(3, -fi5p-p5P, (35u), 

A eq 56 = /T 1 (eV.56> (3) - /^A56> (1) , P VM {2 \ 

eV56 w + pVSQW + 7l V.(D56> (3) - 2d" 1 V.56> (3) ) + 72 V(V.56> (3) ) . (A5) 
Further, the dual, A* q , of A eq is defined by the identity 

(69, A* q F) = (A eq ^, F) V FeP(Q). (A6), 
The formula (2.36) follows immediately from Eqs. (A. 5) and (A. 6). 

Appendix B. Proof of Lemma 3.5. 

By the standard relationship between distributions and the corresponding generalised 
functions, 

©(V/) = - jf dxf(x)^(l%(x) V feV(Q) (B.l) 

and 

Ct ( ?(Vfc) = - jf dxh t (x)^% hJ (x) V ^Pt(O). (B.2) 
On combining these formulae with Eqs. (2.19), (3.21) and (3.22), we see that 

£ dxdx'f(x)f\x')^^E eq (cl%^ = 
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| [s, t]n[s', t'\ V /, f'eviyt) t, s, t', s'eR. (s.3) 



and 



2/3 1 / dxdx'hi{x)hj{x') - (n/i(SijSi m + 6 im Sji-2d 1 S i i5 jm ) + n/2Su6 jm )6(x - x') 

Jq2 OXlOX m 



X 



I [s, t]n[s', t'] I V h, h'eV v (Q), t, s, t', s'eR. (5.4) 

Eqs. (B.3) and (B.4) signify that 



d 



2 



T[ E e*$%(*)C*,.'»&)) ~ 2KS(xx')\[s,t]n[s',t'\] = (5.5) 



dxidx 



and 



<9 



2 



dxidx' m [ Ee ^t^ii( x )(t',s'-,jm( x ')) 

+ 5im5jl — 2d~ L 8 

ilvjm ) + (72<^<5j TO )5(a; - x')x 

Ita.tlnf*',* 7 ]!] =o. (5.6) 



We now invoke our assumptions that the tensor f, and hence Q 3 \ is symmetric, that 
the equilibrium two-point functions of Q 1 ' and are translationally and rotationally 
invariant and that they are of zero range. It then follows that these functions take the 
following forms. 

Ee q {Ci%(x)Ct', 8 '-A x ')) = A (* ~ (B.7) 

and 

Ee*($%(*)&% m &)) = B ( x - + U«) + C(x - x')5 u S jm , (B.8) 

where A, B and C are generalised functions on O that depend on £, s, t' and s' . It follows 
from Eqs. (B.5) and (B.7) that 

A [A(x - x') - 2kS(x - x') | [s, t]n[s', t'} \] = (B.9) 

and from Eqs. (B.6) and (B.8) that 

A[B(x-x') -2p- 1 <y 1 5(x-x')\[s,t]r\[s',t'}\] =0 (BAO) 

and 

A[C{x-x') - 2(3- 1 (- f2 - 2d~ 1 -f 1 )5(x - x')\[s,t]n[s' ,t']\] =0. (5.11) 
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In view of Schwartz's point support theorem [6, Th. 35], these last three equations signify 
that 

A(x-x') = 2K8(x-x')\[s,t]n[s',t'}\, (B.12) 
B(x-x') = 2fi- 1 <y 1 5(x - x')\[s,t]n[s' ,t']\ (S.13) 

and 

C(x - x') = 2/T 1 ( 72 - 2d~ 1 -f 1 )5(x - x') \ [s, t]n[s', t'}\. (£.14) 

Eq. (3.25) now follows from Eqs. (B.7) and (B.12); and Eq. (3.26) follows from Eqs. 
(B.8), (B.13) and (B.14). Finally, in view of the observation following Eq. (3.24), a 
parallel treatment of the other components of the two-point functions of £ reveals that 
they all vanish. 

References 

[I] A. Einstein: Ann. Phys. 11 (1903), 170; 17 (1905), 549 
[2] L. Onsager: Phys. Rev. 37 (1931), 405; 38 (1931), 2265 

[3] L. D. Landau and E. M. Lifschitz: Fluid Mechanics, Pergamon, Oxford, 1984 

[4] G. L. Sewell: Lett. Math. Phys. 68 (2004), 53 

[5] G. L. Sewell: Rev. Math. Phys. 17 (2005), 977 

[6] L. Schwartz: Theorie des distributions, Hermann, Paris, 1998 

[7] G. L. Sewell: J. Phys. A 41 (2008), 382003 

[8] R. F. Streater and A. S. Wightman: PCT, Spin and Statistics, and All That, Benjamin, 
New York, 1964 

[9] G. Grinstein, D. H. Lee and S. Sachdev: Phys. Rev. Lett. 64 (1990), 1927 

[10] J. R. Dorfman, T. R. Kirkpatrick and J. V. Sengers: Ann. Rev. Chem. Phys. 45 
(1994), 213 

[II] H. Spohn: J. Phys. A 16 (1983), 4275 

[12] B. Derrida, J. L. Lebowitz and E. R. Speer: J. Stat. Phys. 107 (2002), 599 

[13] L. Bertini, A. de Sole, D. Gabrielli, G. Jona-Lasinio and C. Landim: J. Stat. Phys. 
107 (2002), 635 

[14] D. Ruelle: J. Stat. Phys. 98 (2000), 57 

[15] L. Accardi, A. Frigerio and J. T. Lewis: Publ. Res. Inst. Math. Sci. 18 (1982), 97 



23 



